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Abstract 

> . 

C<") . Each scheme of state reconstruction comes down to parametrize the state of 

a quantum system by expectation values or probabilities directly measurable in 
• an experiment. It is argued that the time evolution of these quantities provides 

an unambiguous description of the quantal dynamics. This is shown explicitly 
for a single spin s, using a quorum of expectation values which contains no 
redundant information. The quantum mechanical time evolution of the system 
is rephrased in terms of a closed set of linear first-order differential equations 
coupling (2s + l) 2 expectation values. This 'realization' of the dynamical law 
refers neither to the wavefunction of the system nor to its statistical operator. 

q-( For a quantum system with statistical operator p, it is straightforward to determine 

the expectation value (A) of an operator A according to 

><: (A) p = Tr \Ap 

Methods of state reconstruction 0, 0, solve the inverse problem: the unknown 
state p of the quantum system is expressed as a function of the expectation values 
of properly chosen observables Q J , which constitute a quorum Q. The resulting 
equivalences 

P {{Q 3 ),jel}, Q j eQ, (2) 

with j taking on values from a discrete or continuous set I of labels, are more than 
just mathematical beauties — they have been used in the laboratory to reconstruct 
correctly states of various quantum systems || ||, 0]. For example, the state of an ion 
in a Paul trap has been identified || by a method realizing (f2|) on the basis of Wigner 
functions. 

The purpose of the present paper is to point out that a parametrization as in (PJ) 
suggests a conceptually interesting way to describe the time evolution of a quantum 
system without invoking its density matrix or wave function. Instead, only directly 
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measurable quantities, that is, expectation values of hermitean operators are involved. 
The argument will be given in general terms first, specifying neither the system at 
hand nor a particular method of state reconstruction. In the main part of the paper 
the example of a single spin s is worked out explicitly, followed by a discussion putting 
the results into perspective. 

The 'quantum mechanical Liouville' equation 0, 

describes the time evolution of a quantum system with Hamiltonian operator H and 
statistical operator p. The state po at time to is transported to pi at time t\ along 
a smooth path in state space. Suppose that the operators {Q- 7 } provide a quorum 
Q for all possible states of the system at hand. Then, each state pt on the path 
between p and p\ is characterized uniquely by the set of expectation values {(Q j )t}- 
In other words, the path p t in state space has an unambiguous image in the space of 
expectations {(Q- 7 )}. This path is expected to arise as the solution of a dynamical law 
in this space [fLOU : 

j t (&)t = VL({{&) t }) , (4) 

where the function T> will depend on both the Hamiltonian H of the system and the 
quorum Q. Subsequently, time dependent expectation values of arbitrary operators 
A can be expressed in terms of the {(Q-')t} simply by using Eq. (|2|) to eliminate p t 
in favor of the elements of the quorum. 

In the following, an explicit form of Eq. ([|) will be derived for a spin s using 
a nonredundant quorum. Quantum mechanically, the spin is described by a vector 
operator S = hs, the components of which satisfy the commutation relations of the 
algebra sw(2): [s x , s y ] = is z , . . . These operators act irreducibly in a Hilbert space TL S 
of (complex) dimension (2s + 1). The standard basis of the space Ti s is given by the 
eigenvectors of the z component of the spin, S z = n z ■ S, and they are denoted by 

n z ), -s < p < s [|11|. 

Observables are represented by hermitian operators, A* = A, all of which are linear 
combinations of polynomials in the operators s x , s y and s z of degree 2s at most. The 
ensemble of all hermitean operators acting on 7i s can be considered as a vector space 
A s of dimension N s = (2s + 1) 2 . One well-known basis |T2| for these operators is given 



by the multipoles K im , < / < 2s, — / < m < I. They are associated with the group 
SU(2): as a tensorial set, they transform in a simple way under rotations. In terms 
of those, a hermitean operator A can be written as 

7 — ^ \ " A Im 



- A K lm , (5) 



1 Im 



with a unique set of expansion coefficients A 

Given the multipoles K[ mi a second set of N s operators K is known to exist such 
that 



2s 



-Tr 



<W > (6) 
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where the trace taken in the Hilbert space 7i s defines the scalar product of two op- 
erators. The operators {K' m } constitute a second, dual basis of the space A s , being 
determined unambiguously |13[ by the original basis. In the case of multipoles, the 



elements of the dual basis are known explicitly: K lm = K^ m . Using (El), one easily 



determines the coefficients A in the expansion 



A lm = Tr 



AK lm . (7) 



There is a second expansion for selfadjoint operators A in terms of the dual basis. In 
analogy to Eqs. (|5]) and (^) one can write any operator A as a linear combination of 
the operators K lm with coefficients A^ m = Tr [yl^ m ] divided by (2s + 1). 

Here the purpose is to express the dynamical evolution of a quantum spin in 
terms of expectation values. To this end, a basis different from the multipoles 
will be used. Denote the eigenstate of the operator n • s along the direction n = 
(sin t? cos (p, sin $ sin (p, cos and with eigenvalue s by 

|n) = exp[— i m((p) ■ s] \s, n 2 ) , (8) 

where 111(93) = (— sin <p, cos (p, 0): the state |n) is obtained from rotating the state 
\s,n z ) about the axis m(<£>) in the xy plane by an angle In this way, a coherent 



state |n) is associated to each point of the surface of the unit sphere |fL4|l . The ensemble 
of all coherent states provides an overcomplete basis of the Hilbert space 7i s . 

The density matrix p of a spin s is determined unambiguously if one performs 
appropriate measurements with a traditional Stern-Gerlach apparatus. Distribute 
N s = (2s + l) 2 axes n^, — s < /i, v < s, over (2s + 1) cones about the z axis with 
different opening angles such that the set of the (2s + 1) directions on each cone is 
invariant under a rotation about z by an angle 27r/(2s + 1). As shown in an 
unnormalized statistical operator p is then fixed by measuring the (2s + l) 2 relative 
frequencies p s (n MI/ ) = (n MI/ |p|n MI/ ), that is, by the expectation values of the statistical 
operator p in the coherent states |n MJ ,). In other words, the projection operators 

QtMu = |n MI/ )(n^| , s<fi,v<s, (9) 

constitute indeed a quorum Q for a spin s. This fact, when used for state reconstruc- 
tion, defines an optimal method since exactly (2s + 1) 2 numbers have to be determined 
experimentally which equals the number of free (real) parameters of the (unnormal- 
ized) hermitean density matrix p. For the following, it is convenient to replace the 
labels p, and v by a single index n = (p,u), say, with 1 < n < N s . 

The set of all unnormalized hermitean density matrices for a spin s is just the set 
of all hermitean operators acting on the Hilbert space 7i s . Therefore, the quorum Q 
automatically provides a basis for hermitean operators A: 



A 



I N s 



2* + * n =l 



-J2A n Q n , A n = Tr L^LQ^J , (10) 

-L 1 



and the expansion coefficients A n involve operators Q n dual to the elements of the 
original basis. In analogy to the example involving multipoles, there are orthogonality 
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relations, 



1 



Tr 



5 n 



1< n, n' < N« 



2s + 1 

and a second expansion for hermitean operators is available: 



111 



A 



1 



2s + 



Tr 



A Qn 



71=1 



121 



Let us now consider the properties of the statistical operator p when expanded in 
the basis Q n dual to the original quorum, 



2s 



(13) 



71=1 



where the coefficients P n = Tr pQ n = (n n |p|n n ) satisfy 



< P„ < 1 



1 < 71 < iV, 



(14) 



Each of the N s numbers P n has a value less or equal to one due to the normalization 
of the density matrix, Tr [p] = 1, and since p is a positive operator, the P n are 
non-negative throughout. This is a unique and essential feature of the basis {Q n }— 
the expansion coefficients of p with respect to neither the original basis {Q n } nor 
the multipole bases {Ki m } or {K' m } have this property. The interpretation of the 
coefficients P n — to measure the value s along the axis n ra — is clearly compatible with 
fll4). It is important to note that, although each of the P n is a probability, they do 



not sum up to unity: 



< J2 R n < ( 2S + !)' 



(15) 



71=1 



This is due to the fact that they all refer to different orientations of the Stern-Gerlach 
apparatus, being thus associated with the measurement of incompatible observables, 



Q n , Q n > ± , 1 < n, ri < N s 



(16) 



since the scalar product (n n |n n /) of two coherent states is different from zero. The sum 
in ( 15) cannot take the value (2s + l) 2 since this would require a common eigenstate 



of all the operators Q n which does not exist due to (|16|). By an appropriate choice 
of the directions n ra (all in the neighborhood of one single direction no, say), the sum 
can be arbitrarily close to (2s + l) 2 for states 'peaked' about n . Similarly, the sum 
of all P n cannot take on the value zero since this would require a vanishing density 
matrix which is impossible |T5 |. If, however, considered as a sum of expectation values, 
there is no need for the numbers P n to sum up to unity. Nevertheless, they are not 
completely independent when arising from a statistical operator: its normalization 
implies that 

1 



N s 



Tr 



2s + 



-Y Tr 
l±i 



Q' 



1 



(17) 



71=1 
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which turns one of the probabilities into a function of the (2s + l) 2 — 1 = 4s(s + 1) 
others, leaving us with the correct number of free real parameters needed to specify 
a density matrix. 

It is useful to visualize the description of a density matrix by the numbers P n in ge- 
ometrical terms. Consider the linear space A s of dimension N s , each axis being associ- 
ated with one projector Q™ and the coefficient P n . Since Tr [Q n Q n i] = |(n n |n n /)| 2 ^ 0, 
this is not an orthonormal basis of A s , and neither is its dual {Q n }. According to 
(|13D and (|l4l) each statistical operator determines a point P with components P n in an 
iV,,- dimensional parallelepiped [jT5| . Eq. ( fPTD may be understood as a scalar product 
of P and the vector E with components E n = Tr [i?Q n ] = Tr [Q n ] , where E denotes 
the unit operator in A s . Thus, the points P which correspond to normalized density 
matrices are necessarily located in the intersection 1Z of the parallelepiped with a 
hyperplane of dimension 4s(s + 1) in A s . However, not all points in 1Z are associated 
with a density matrix. To see this, imagine the quantum system to be in an eigenstate 
|s,n no ) of the projection operator Q no , say. Then, the corresponding probability P no 
has the value one, and all the others are smaller than one. This is the only point of 
the unit cube with P no = 1 associated with a density matrix but one constructs easily 



other points satisfying ([17|) . 

Let us turn to the dynamics of the quantum system expressed by the probabilities 
P n . Their time derivative, dP n /dt = (n n \dp/dt\n n ) , is determined unambiguously by 
Eq. (^). Using the expansion flT5| ) it is easy to express the resulting equations in the 
form (H). A closed set of equations for the variables P„(t) follows from plugging (|T3| ) 
into the right-hand-side of Eq. (§) and taking the expectation value in the state |n n ): 



Ns 



4p»(*) = lT,^n\[Q n ',H]\n n )PAt) 

n'=l 



Thus, the spin dynamics has been expressed entirely in terms of the N s variables P n : 
this equation, the explicit form of (|j) for a single spin s, provides the main result of 
this paper. In fact, the dynamics is consistent with (|T7|): multiply (|T8|) by Tr [Q n ] 
and sum over all values of n: 



n=l 



d 



Tr [Q«j-P n = £ Tr \Q n ',H 



Tr 







(19) 



using ( TTS] ) and expanding the identity as E = J2 n Tr [Q n ]Q n . Consequently, the time 
evolution of the quantum system can be represented by a point moving in the domain 
71, with a trajectory determined by (0). Eq. ([18]) will be called the 'expectation- value 
representation' of the equation of motion (0). 

Let us point out some properties of the time evolution of the spin s when given 
in the expectation-value representation. The dynamical law (|18|) is a closed set of 
linear equations for the N s real variables P n : the time derivatives dP n /dt at time t 
are expressed entirely in terms of the probabilities P n at that time. Introduce a real 
A", x iVo matrix M with entries 



M n "' = ^(n n |[Q"',F]|n n ) = ^Tr 



H[Q n ,Q n ' 



(20) 
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using the cyclic property of the trace. Then, one can rewrite the dynamics Jig) and 
the constraint ( |T7| ) as 

^ = MP(i), P(t )-E = l. (21) 

Therefore, the quantum dynamics of a spin s is equivalent to that of a classical 
dynamical system with N s degrees of freedom, constrained to move in a certain region 
1Z to be considered as its phase space. For an eigenstate of the Hamiltonian H with 
eigenvalue and density matrix pS k > = |efc)(efc|, one has [p^ k \ H] = 0; hence, the 
flow generated by M in 1Z has precisely (2s + 1) 'fixed points' with coordinates P^. 
A more detailed study of the flow in ( pT| ) will exploit the existence of a metric in 
the space A s , induced by the transformation between the original basis and its dual: 
Q n = (2s + l) -1 J2 n ' Gnn'Q™ • This metric, G nn > = Tr [Q n Q n '}, can be shown to define 
a positive definite quadratic form. 

Somewhat surprisingly, the function T> introduced in (^) is linear in the variables 
P n which, in turn, are linear functions of the density matrix p. Therefore, the convexity 
of the state space, p (A) = (1 - A)p (a) + Ap (6) , < A < 1, turns into: 

p(A) = p(a) + A fp{b) _ p(a)\ ^ ( 22 ) 



tracing out a straight line in the space of expectations. 

In the expectation-value representation, the distinction between different pictures 
of quantum mechanics (Schrodinger, interaction, Heisenberg) no longer applies since 
the quantum dynamics has been expressed entirely in terms of observable quantities. 
Furthermore, in this representation there is no temptation to attach a wave function to 
an individual quantum spin as is done in the 'individual' interpretation of quantum 
mechanics ||17|| . From the outset, the involved probabilities make sense only when 
referring to an (infinite) ensemble of spin systems prepared identically. Therefore, a 
vector P is associated rather with the procedure of state preparation than with an 
individual microscopic system. Turned around, the expectation-value representation 



seems to favour the 'statistical interpretation' of quantum mechanics [17]. 

Conceptually, the 'realization' introduced here differs from other formulations of 
quantum mechanics 'without wave function' such as the phase-space representation 
through Wigner functions, be it for a particle JOJ or a spin |TJ|. The occurrence of 
negative values is characteristic of 'quasi-probability' distributions a la Wigner, ex- 
pressing the impossibility that position and momentum simultaneously have definite 
values. The expectation-value representation allows one to rephrase the quantal dy- 
namics entirely in terms of directly observable and non-negative quantities, defined on 
(2s + l) 2 points of the sphere. This construction can be thought of as discretizing the 
phase space of the classical spin and associating probabilities (which, however, do not 
sum up to unity) to the individual points. There is a conceptual link to the 'prob- 
ability representation' for both quantum particles and spins fl20| , f2~l|l which is based 
on positive smooth distributions on the classical phase space of the underlying sys- 
tem. It provides, however, a highly redundant description, while the expectation- value 
representation works with nonredundant information only. 
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From a general point of view, Eqs. ([M]) and ( O) provide the basis for a sym- 
bolic calculus comparable to the Wigner formalism |2"3| or to the coherent state 
representation |24| of quantum mechanics. Once the vector P(to) associated with a 



quantum state p(io) is known, one can extract the time evolution of arbitrary observ- 
ables A(t) from P(t) without ever invoking p{t). The details of this calculus based on 
the expectation-value representation will be developed elsewhere. Furthermore, it is 
not difficult to generalize the present approach to non- autonomous quantum systems 
described by explicitly time dependent Hamiltonian operators H(t). 

In sum, the expectation-value representation of quantum mechanics, as derived 
from (|K|), is equivalent to any other representation. The statistical operator p of a 
quantum spin s is represented by a point on a manifold in an ^-dimensional space 
parameterized by probabilities or, equivalently, expectation values. Its time evolution, 
Eqs. ( |T8"D or (pip, is governed by a linear and autonomous classical flow. Therefore, one 
can describe the quantum dynamics as a smooth trajectory in the space of expectation 
values. With all unobservable elements eliminated from the theory, the expectation- 
value representation provides an appealing explicit realization of Schrodinger's remark 
to consider the wave function as a "Katalog der Erwartung." 
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